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ABSTRACT 



The blazar 1ES 0229+200 is a high frequency peaked BL Lac object with a hard TeV spectrum extending to 10 TeV. Its unusual 
spectral characteristics make it a frequently used probe for intergalactic radiation and magnetic fields. With new, simultaneous ob- 
servations in the optical, ultraviolet (UV) and X-rays, the synchrotron emission is probed in great detail. The X-ray emission varies 
by a factor of « 2 in 2009, while being rather stable in 2010. The X-ray spectrum is very hard (T ss 1.8) and it shows an indication 
of excess absorption above the Galactic value. The X-ray emission is detected up to ~ 100 keV without any significant cut-off, thus 
1ES 0229+200 belongs to the class of extreme blazars. The simultaneous measured, host galaxy- and extinction-corrected optical and 
UV fluxes illustrate that the cut-off of the low energy part of the synchrotron emission is located in the UV regime. The minimum 
energy of the electron distribution has to be rather high to account for this cut-off. This implies that there is a narrow-band energy 
distribution function of radiating electrons, which is responsible for the unusually hard TeV spectrum. Other extreme blazars have 
similar synchrotron peak frequencies but much softer TeV spectra, hence 1ES 0229+200 has one of the highest inverse Compton (IC) 
peak frequency and the narrowest electron distribution among the extreme blazars known to date. 
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1. Introduction 

The high-frequency peaked BL Lac object 1ES 0229+200 is 
located at gnnon = 2 h 32 m 48.62 s , r) J2 ooo = +20 ° 17' 17.45" 
dRector et al.ll2()()3l) and has a redshift of z = 0.14 dWoo et al.l 
120051) . 

High frequency peaked BL Lac objects (HBL) are charac- 
terized by two peaks in their spectral energy distribution (SED) 
which are located in the UV-X-ray and the GeV-TeV band, re- 
spectively. T hese are commonly inte rpreted in terms of leptonic 
models (e.g. lMarscher & Gearll 19851) as synchrotron and inverse 
Compton (IC) emission from a population of relativistic elec- 
trons upscattering their self-produced synchrotron photons (syn- 
chrotron self-Compton (SSC) models). 

1ES 0229+200 was d iscovered in the Einstein IPC Slew 
Survey (Elvi s et al .1 119921) . and classified as a high-frequency 
peaked BL Lac obj ect based on its X-ray to radio flux ratio 
dGiommi etaD 19951) . 

The VLA observations of 1ES 0229+200 reveal a core flux 
of 51.8 mJy and show curved jets with an extension of ~ 30" at 
1.4 GHz and jet po sition angles of P.A.= -10° and P.A.= 180° 
dRector et alj|2003l) . 

1ES 0229+200 is not detected in the high energy y-ray range 
(100 MeV < E < 100 GeV) by Fermi/LAT in two years of obser- 
vations and hence it is not mentioned in the second Fermi catalog 
dAbdo et al.ll20TTh . 

In 1996, it was originally pr edicted to be a pote ntial VHE 
y-ray source based on its SED (iStecker et al.l [l996). however, 
Whipple, HE GRA, a nd Milagro have o nly reported upp er hmits 
dHoran et all d2004l) . lAharonian etail (2004). Wi lliams! d2005l) '). 
Very high-energy (VHE, E>100 GeV) emission up to 10 TeV 
was first detected with the High En ergy Stereoscopic System 
(H.E.S.S.) in 2006 (lAharonian et al.ll2007l) . In this study, a hard 



spectrum with a photon index of F - 2.5 + 0.19 stat + 0. 1 5V5 was 
reported. Besides 1ES 1426+428 dAharonian et al.l .20031) ). 1ES 
0229+200 is the only source at redshift z > 0.1 whose spec- 
trum was measured up to this high energy. 1ES 0229+200 has a 
very hard int rinsic VHE spectrum, he n ce it is well-suite d to EB L 
studies (e.g. lAharonian et al. (2007), Kneiske & Dold ((20ll). 
Its spectral characteri stics have been used to p r obe intergalactic 
magne tic fields (e.g. iNeronov & Vovkl d2010h . iTavecchio et al.l 
(2010)). At the same time, these hard spectra are challenging 
for blazar models. We have launched dedicated, simultaneous 
multi-wavelength observations with XMM-Newton (X-ray, UV) 
and ATOM (optical) to determine the broad-band spectra in more 
detail. 



2. Multi-wavelength observations and data analysis 

2.1. X-ray data from XMM-Newton, Swift, and RXTE 

XMM-Newton observations of 1ES 0229+200 were carried out 
on August 21 and 23, 2009 for 23 and 28 ks, respectively. The 
observations were conducted with MOS 1 and PN in full imag- 
ing mode and MOS2 in timing mode, all with a thin filter. The 
two grating spectrometers onboard XMM-Newton RGS 1, 2 were 
also used to acquire. The data analysis of the XMM-Newton ob- 
servations was performed with SAS v.9.0. The data from the 
X-ray instruments were reprocessed as described in the SAS 
user-guidd_l Both observations are influenced by short soft pro- 
ton flares at the end of each observation, detected in the high 
> 12 keV emission for the PN (> 10 keV for MOS) detector. 
Therefore, the good time intervals were determined using a cut 



1 http://xmm.esac.esa.int/external/xmm_user_support/documentation 
/sas_usg/USG/ 
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Instrument 


Time 


ObsID 


Energy range 


Photon index 


^2-10keV,deab S (erg Cm" 2 S _1 ) 


XMM-Newton 
Swift /XRT 
Swift /XRT 
RXTE /PC A 
BeppoSAX 
ROSAT 
Einstein 


21., 23 August 2009 
5.,7.,8. August 2008 
19.0ct. - 23. Nov. 2009 
Uan. - 13. Oct. 2010 
16. July 2001 


604210201,604210301 
31249001-31249003 
31249004-31249019 
95387 
51472001 


0.1-15 keV 
0.2-10 keV 
0.2-10 keV 
3-60 keV 
0.1-50 keV 
0.1-2.4 keV 
0.4-4.0 keV 


1.84 ± 0.02 
1.87 ±0.05 
1.73 ±0.03 
1.92 ±0.05 
1.99 ±0.05 


(9.0 ± 0.1) x 10~ 12 
(1.03 ±0.04) x 10 -11 
(1.47 ±0.04) x 10-" 
(1.23 ±0.04) x 10"" 
1.5 x lO -11 
(4.5 ± 0.6) x 10~ 12 at IkeV 
(7.6 ± 2.2) x 10~ 12 at IkeV 



Table 1. X-ray observations of 1ES 0229+200 showing the results for the annual binning of X MM-Newton, Swi ft, and RXTE, as wel l 
as historical observations from Einstein (lElvis et al.lll992l) . ROSAT (Brinkm ann et al.lll995l) and BeppoSAX (Don ato et al.l 12005). 
The measured fluxes are given for the energy range between 2 and 10 keV. For the ROSAT and Einstein observations, the fluxes are 
determined at an energy of IkeV. The absorption derived from the high precision XMM-Newton spectra of Nh — 1.1 x 10 21 cm 2 was 
used for the XMM-Newton, Swift, a nd RXTE spectra to determine the unabsorbed flux. The absorption detected in the BeppoSAX 
spectrum was JV H = 1.0 x 10 21 cm 2 dDonato et alj|2005l) . 



at 0.4 cts/s for PN (0.35 cts/s for MOS) for the high-energy count 
rate resulting in an effective exposure of 7.3 ks for PN (23.1 ks 
for MOS) for the first and 12.1 ks for PN (24.4 ks for MOS) for 
the second pointing. No significant variations were found in ei- 
ther observation using different binning down to a timescale of 
100 (10) seconds for the imaging mode (timing mode) during 
these pointings. 

For the imaging mode, spectra were extracted from a source 
region of 30" radius for the PN (80" for MOS) around the po- 
sition of 1ES 0229+200 and background regions with the same 
radii, on the same chip of the detectors that contained the source 
region. No significant pileup was found in the spectra of the 
source regions using the tool epatplot. The RMF (redistribu- 
tion matrix file) and ARF (ancillary response file) are calculated 
for each spectrum using the tools rmfgen and arfgen. 

The grating spectrometers RGS 1, 2 onboard XMM-Newton 
measure in the energy range 0.35 to 2.5 keV. The data are anal- 
ysed using the tool rgsproc. 

From October 19 to November 23, 2009, 16 Swift obser- 
vations were targeting on 1ES 0229+200. During these point- 
ings with observatio ns of 0.5 - 4 ks each, the XRT detector 
(Burr ows et al J 120051) was operated in windowed-timing (WT) 
mode in the energy range 0.2 - 10 keV. In 2008, 1ES 0229+200 
was observed on August 5, 7, and 8 in photon-counting (PC) and 
WT mode. However, only the PC mode observations are taken 
into account in the analysis, since only 120 s were observed in 
WT mode resulting in an insufficient amount of counts. For the 
Swift spectral analysis, XRT exposure maps were generated with 
the xrtpipeline to account for some bad CCD columns that 
are masked out on-board. The masked hot columns appeared 
when the XRT CCD was hit by a micrometeoroid. Spectra of 
the Swift data in PC-mode have been extracted with xselect 
from an annulus region with a radius of 0.8' at the position of 
1ES 0229+200, which contains 90% of the PSF at 1.5 keV. The 
background was extracted from a circular region with radius of 
3' near the source. For the WT-mode, boxes (~ 1.6' x 0.3') cov- 
ering the region with source photons and a background region 
of similar size were used to extract the spectra. The auxiliary 
response files were created with xrtmkarf and the response 
matrices were taken from the Swift package of the calibration 
database caldb 4 . 1 . ]Q. 

X-ray observations with t he Proportional Counter Array 
(PC A) detector onboard RXTE dBradt et all 19931) were obtained 
in the energy range 2-60 keV from January 1 to October 13 
2010 with exposures of 1 - 2 ks per pointing. Only RXTE /PCA 
data of PCU2 and the top layer 1 were considered to ensure the 

2 http://heasarc.gsfc.nasa.gov/docs/heasarc/caldb/caldb_intro.html 



highest signal-to-noise ratio. The data were filtered to account 
for the influence of the South Atlantic Anomaly, tracking off- 
sets, and electron contamination using the standard criteria rec- 
ommended by the RXTE Guest Observer Facility (GOF). For 
the count rate of ~ 1 cts/s for this observations, the faint back- 
ground model, provided by the RXTE GOF was used to gener- 
ate the background spectrum with pcabackest and the response 
matrices were created with pcarsp. 

The second instrument onboard RXTE, the HEXTE (High 
Energy X-ray Timing Experiment) onboard RXTE takes data in 
the energy range 15 to 250 keV. Since 2006, the HEXTE cluster 
A has been operating only in ON-source mode and since the end 
of 2009 the cluster B is permanently operating in OFF-source 
mode. For the spectral analysis, the cluster B data are used as 
background information for the cluster A data. The sum of all 
observations from January 1 until October 13, 2010 show no 
significant signal from 1ES 0229+200 in this energy range. 



2.2. Additional X-ray sources in the XMM-Newton field of 
view 

The source detection tool edetect_chain revealed 20 point 
sources in the field of view of the PN detector and one source 
that is extended beyond the PSF of the instrument. Two point 
sources are coincident with IRAS sources, four point sources 
are coincident with NVSS s ources, and 1 5 with stars from the 
Guide Star Catalogue (GSC, lLasker et ail d2008l) ) (see table in 
online supplement). The remaining point sources do not have a 
counterpart in these catalogues and remain unidentified. 

The extended source XMMU 023318.0+201237 is located 
at ttfiooo = 2 h 33 m 18.05 s , 5 J2 ooo = +20°12'37.19" and is very 
close (58" offset) and likely connected to a point-like source 
XMMU 023315.5+201323 that is positional coincident with an 
infrared source IRAS 02304+2000 with coordinates anooo = 
2 h 33 m 14.6 s , fe non = +20°13'30" and a flux of 0.4 Jy at 60 pm 
dMoshirll9 90). This infrared (IR) source has a radio counterpart 
NVSS023314+201330 that has an elliptical shap e in the NVSS 
sky m ap and a flux of 12.4 ± 1.1 Jy at 1.4 GHz dCondon et al.l 
1998). One of the R-band observations by ATOM cover the re- 
gion of this source and an object with 16.5 mag is measured. 
A closer look yields the detection of two very close (4" separa- 
tion) point sources of which the brighter one can be identified 
with USNO-B1.0 1102-0028956. 
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2.3. UVdata from XMM-Newton/OM and Swift/UVOT 

The optical monitor (OM) dMason et all 1200 lb onboard XMM- 
Newton observed 1ES 0229+200 in the filters UVM2 (231 nm), 
UVW1 (291 nm), and U (344 nm) simultaneously with the X- 
ray telescope. The analysis of these data was performed with the 

analysis described in the webpag^E 

The UVOT instrument (Roming et al. 2005) onboard Swift 
measures the UV and optical emission in the bands UVW2 (188 
nm), UVM2 (217 nm), UVW1 (251 nm), U (345 nm), B (439 
nm), and V (544 nm) simultaneously with the X-ray telescope 
with an exposure of 0.4-3 ks ea ch. The instrumen tal magnitudes 
and the corresponding flux (see lPoole et al .1120081 for the conver- 
sion factors) are calculated with uvotmaghi st taking into ac- 
count all photons from a circular region of radius 5" (standard 
aperture for all filters). An appropriate background was deter- 
mined from a circular region of radius 5" near the source region 
without contamination of sources. 

2.4. Optical data from ATOM 

The 75-cm telescope ATOM (lHauser et al.ll2004l) . located at the 
H.E.S.S. site in Namibia, monitored the flux in the different fil- 
ters B (440nm), V(5 50 nm), R (640 nm), and I (790 nm) accord- 
ing to lBesselil ( 19901) . The obtained data were analysed using an 
aperture of 4" radius. Photometric calibrat ion was done using 
the standard fields S A 1 1 3 and S A 95 from|Landolt| dl992l) . 

3. Spectral data in the synchrotron range 

All XMM-Newton, Swift, and RXTE spectra are binned with at 
least 25 counts and xspec vl2 . 5 was used for the spectral anal- 
ysis. For the XMM-Newton spectra, the energy ranges were re- 
stricted to 0.1 - 10 keV for MOS and 0.2-15 keV for PN 
following the suggestion of the calibration informatioifl The 
MOS1 and PN spectra are fitted simultaneously using a con- 
stant parameter to account for the different normalizations. A 
power-law model of the form F(E) = Nq{E I Eq)~ y was used to 
fit the X-ray spectra taking into accoun t the Galactic absor ption 
ofNu = 7.9 x 10 20 cm 2 (LAB Survev. lKaiberla et all (120051) 1. As 
can be seen in Fig.QJ, the residuals of the fit to the X-ray spectra 
deviate from the expected form, which implies that their is either 
additional absorption or different spectral characteristics. 

Fitting a power-law model with unconstrained absorption re- 
sults in a much better description of the data with % 2 /dof = 
979.7/903 = 1.08. The fit parameter and the goodness of each 
model is given in table|2] Residuals are shown in panel b of Fig. 

m 

Alternatively, the deviations shown in Fig. [T] could be 
avoided by taking into account the Galactic absoption and fit- 
ting a broken power-law. This fit results in photon indices of 
Ti = 0.8±0.3,F 2 = 1.78 ±0.01, and a break at 0.58 ± 0.09 keV 
(X 1 /dof = 935/887). However, the low energy extrapolation of 
the X-ray spectrum would not fit the UV-optical range, hence 
this model is disfavoured. 

To identify the location of the additional absorption, we 
test the hypotheses that either the total A^h is located at red- 
shift z=0 or that the additional absorption is located at redshift 
z=0.14. The resulting models with two different redshifts are in- 
distinguishable and have similar goodnesses of fit (at z=0.14: 



X 2 /dof = 961.3/888; at z=0: x 2 /dof = 959.3/888). Hence the 
location of the additional absorption could be either intrinsic to 
1ES 0229+200 or in the line of sight to the observer or in the 
Milky Way. We note that both a local enhancement of Galactic 
N H by 25% along the line of sight to 1ES 0229+200, as well 
as an intrinsic column within the source of 2.9 x 10 20 cm 2 , are 
plausible. 
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3 http://xmm.esac.esa.int/sas/current/documentation/ 
threads/omLstepbystep.shtml 

4 XMM-SOC-CAL-TN-0018: 
http://xmm.vilspa.esa.es/docs/documents/CAL-TN-0018.pdf 



Energy (keV) 

Fig.l. XMM-Newton MOS1 (black) and PN (red) spectrum of 
1ES 0229+200 from August 21, 2009. The spectra can be well 
fit with a power-law model taking into account an absorption 
larger than the Galactic absorption. In panel a, we plot the resid- 
uals for a power law considering the Galactic absorption as fixed 
parameter, and in panel b residuals for a power law with a free 
absorption. 



The XMM-Newton/RGS spectra show no significant line 
emission and the continuum spectra are well-described by a 
power law with additional absorption comparable to the PN and 
MOS spectra. 

The spectrum of the extended source detected in the PN and 
MOS observations close to 1ES 0229+200, was extracted from a 
region of radius 1 .7' for the PN detector ( 1 ' for the MOS detector 
due to the CCD gaps) and reveals 90% and 30%a faint source 
with a flux of » 1 x 10~ 13 erg crrT 2 s . We were able to fit 
the spectrum with a binning of at least 15 photons per bin with 
a fixed Galactic absorption and a power law model with T = 
2.0 ± 0.2 (x 2 /dof = 73/96) and slightly better by a thermal 
model (mekal) with kT -4 + 2 keV (x 2 /dof = 87/96). 
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Date 


N a 
(cm 2 ) 


r 


Norm. 

(ph cm" 2 s -1 keV" 1 ) 


Norm, factor 
(PN) 


X 2 /dof 


F2—W keV 

(erg cirT 2 s -1 ) 


21. August 2009 


7.9 x 10 2U 
(1.08 + 0.04) x 10 21 


1.72 + 0.01 
1.84 + 0.02 


(2.44 ± 0.03) x 10- J 
(2.73 ± 0.05) x 10~ 3 


0.990 
1.004 


1125.5/889 
959.3/888 


(9.6 + 0.1) x 10~ 12 
(8.9 + 0.1) x 10~ 12 


23. August 2009 


7.9 x 10 2U 
(1.06 ± 0.03) x 10 21 


1.72 + 0.01 
1.83 + 0.02 


(2.49 ± 0.02) x 10~ J 
(2.77 ± 0.04) x 10~ 3 


0.974 
0.989 


1376.9/1064 
1133.7/1063 


(9.8 + 0.1) x 10 12 
(9.1+0.1) x 10~ 12 



Table 2. XMM-Newton fit parameter, goodness of fit, and unabsorbed flux resulting fr om a simultaneous fit of the MOS 1 and PN 
X-ray spectra. The absorption was fixed to the Galactic absorption of the LAB survey dKalberla et al.ll2005l) in the first model and 
let free to vary for the second resulting in additional amount of absorption. 



The annual averages of the Swift spectra in 2008 and 2009 
are shown in tableQ] The consecutive pointings of the 2009 data 
obtained between October and November were binned in three 
ten-day intervals each to increase the statistics and are shown 
in Fig. [2] Initially, a simple power law and a free absorption 
was used to fit these spectra and result in the best description 
of the spectral shape. No significant change in absorption was 
detected with values of N n = (1.6 ± 0.4) x 10 21 cm 2 in 2008 
and N H = (1.3 ± 0.5) x 10 21 cm 2 , N H = (0.9 ± 0.5) x 10 21 cm 2 
and N H = (1.3 ± 0.2) x 10 21 cm 2 for the spectra of 2009 binned 
in ten-day intervals. These values are comparable to the one ob- 
tained from the XMM-Newton spectra, which provide the most 
precise determination of the additional absorption. A total col- 
umn of Nh = 1.1 x 10 21 cm 2 was assumed to obtain the pho- 
ton indices and the fluxes (see table [TJ and Fig. [2J. The Swift 
data from August 2008 were analysed and are shown in the left 
panel in Fig. [I] (2454685 JD). I t should be noted here that these 
data were already presented in iTavecchio e t al. (2009). Our re- 
analysis reveals that the integrated flux between 2 and 10 keV is 
lower by a factor of ~ 7, but confirms the spectral indices. This 
is independent of the choice of a higher value of Nh- 

The RXTE spectra have been summed covering about 30 
days to achieve higher statistics. The energy range from 2-3 
keV and 20-30 keV were excluded from the spectral analysis, 
because of instrumental features (e.g. spikes in the background 
files, see websit^t). A simple power-law was used to fit these 
spectra. The resulting photon indices and fluxes can be seen in 
Fig. El 

The historical spectrum of BeppoSAXwas fitted with a power 
law of photon index F = 1 .99 + 0.05 and a free absorption of 
Nh = (10 ± 5) x 10 20 cirT 2 in the energy range 0.1 - 50 keV 
without any indica tion for a cut-off at the high energy end 
(iDonato et alJl2005l) . 

The Swift/BAT sp ectrum taken from the 58 month catalog 
(Baumgar tnerl d2010h ) is well fit by a power law with photon 
index of T = 2.1 ± 0.3 0f 2 ed = 1.02) and result in a flux between 
14 and 195 keV of F = (28.2 + 5 5) x 10" 12 erg cnr 2 s _1 . The 
BeppoSAX PDS spectrum of 2001 (IDonato et al.l2005l) is in good 
agreement with this high energy spectrum of SwiftfBKT, 
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n 
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2454685 
time / JD 



2455100 2455200 2455300 

time / JD 



Fig. 2. The X-ray flux of 1ES 0229+200 varied by a factor of ~ 2 
in October 2009 as shown in the combined lightcurve taken with 
the instruments XMM-Newton (red cirlces), Swift (black open 
squares), and RXTE (blue crosses) in the energy range 2-10 keV. 
For comparison, the Swift observation of 2008 is also shown. 
The lower panel displays the spectral indices derived by fitting 
power-law models to the spectra. 
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Fig. 3. Ultraviolet fluxes measured by XMM-Newton and Swift 
for 1ES 0229+200 in 2008 (average over 3 days) and 2009 
(U: black circles, UVWI: blue open squares, UVW2: red 
crosses, UVM2: brown open diamonds). The flux points at JD 
2455064.6, 2455066.6, 2455 1 17.8, and 2455 1 1 8.9 are shifted by 
+0.4 days for better visibility. Marginal variations are detected 
in the UVWI and UVM2 bands. 



4. Temporal study 

Regular monitoring observations with ATOM starting in 2006 
with 7 observations in 2006, 47 in 2007, 86 in 2008, 61 in 2009, 
and 24 in 2010 show no significant variations (< 1%) in the R 
band over five years of observations. The average brightness is 
m R = 16.39 ± 0.01 mag and m B = 18.38 + 0.02 mag. 

The UV results from XMM-Newton and Swift of 2008 and 
2009 are shown in Fig. [3] The U and UVW2 band emission was 



3 http://www.universe.nasa.gOv/xrays/programs//JATii /pca/doc/ 
bkg/bkg-2009-spikes/ 
6 http://heasarc.gsfc.nasa.gov/docs/SvW/i f /results/bs58mon/ 



stable and does not show significant variation (p x i = 70% and 
p x 2 = 98%, respectively.). The data in the UVWI and UVM2 
bands display a marginally significant drop by ~ 30% and ~ 
20%, respectively, at JD 2455 120, around the time of an increase 
in the X-ray flux (see below). 

The X-ray flux of 1ES 0229+200 instead shows an increase 
by a factor of around two starting at JD 2455120, as can be 
seen in Fig. [2] For all of 2010, only a marginal variation could 
be detected. During 2008-2010, marginally significant spectral 
changes could be detected as variations in the photon index 
(p x i — 1% for a constant) that were uncorrelated with the flux 
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variation (see Fig. |2j. During all epochs, the source can be well- 
described with a power-law model and an absorption in excess 
of the Galactic absorption, as described in section 3. At ener- 
gies above 14keV, measured by Swift/BAT, marginally signifi- 
cant variation is detected in the monthly lightcurve with a prob- 
ability for the fit of a constant of p y i « 1% over the 58 months 
of observatiorFI ( Baumga rtnerl (120 1 Oh ). However, since a long in- 
tegration time is needed, a behaviour similar to that in the 2-10 
keV range would be difficult to detect. 

5. Spectral energy distribution 




^ io- 3 0.01 0.1 1 10 



energy [keV] 

Fig. 4. Synchrotron emission of 1ES 0229+200 with simulta- 
neously measured optical, UV, and X-ray emission by ATOM 
and XMM-Newton of August 21, 2009 (black dots). The opti- 
cal and UV emission is corrected for both the host galaxy and 
Galactic extinction, and the X-ray emission observed by XMM- 
Newton/MOS is corrected for the absorption. The bars in the 
UVW1 and UVM2 bands are explained in section 5.1. The grey 
open squares represent the non-simultaneous Swift/XRT spec- 
trum (corrected for detected absorption) with the highest flux in 
2009. The Swift/BAT 58 month (Dec. 2004 - Oct. 2009) spec- 
trum is shown in the energy band > 10 keV (black crosses). 
The dashed line (butterfly) re presents the BeppoSAX spectrum 
from 2001 (iDonato et al.ll2005l) . The grey upper limits represent 
historical UV observations with GALEX that are extinction cor- 
rected and of origin discussed in the text. 

The simultaneous observations in the optical by ATOM, the 
UV, and X-rays by XMM-Newton are considered to study the 
synchrotron spectrum (see Fig. [4] simultaneous data of August, 
21 2009). 

5.1. Host galaxy of 1ES 0229+200 

The host galaxy of 1ES 0229+200 is an elliptical galaxy with 
a brightness of mh os t,R = 15 .85 ± 0-01 mag and a half-light 
radius of r e = 3.25 + 0.07" (lUrrv et al.1120001) . Other observa- 
tions in the Bessel R-band with th e Nordic Optical Telescope 
(NOT) ( Falom o & Kotilaiiienll 19991) show results with mh os t,R = 
15.76 mag and r e R = 4". The host galaxy profile of 1ES 

7 available on http://heasarc.gsfc.nasa.gov/docs/Swift/results/bs58mo] 




frequency [Hz] 



Fig. 5. Spectral energy distribution of 1ES 0229+200 with 
simultaneous measured optical, UV, and X-ray fluxes, all 
corrected for host galaxy emission, Galactic extinction, and 
Galactic absorption is shown as black data points. The 58 months 
SwiftfBAT spectrum is shown > 10 keV (black crosses). In grey 
(filled and open circles), historical radio and UV data are shown 
and their origin is discussed in the text. The grey uppe r limits 
in the GeV energy range are taken from iDermer et al .] (120101) 
and represent the upper limit in the energy bins 0.1-1 GeV, 
1-10 GeV, and 10 - 100 GeV by Fermi observations (Aug. 
2008 to Sep. 2010). The VHE y -ray spectrum measure d by 
H.E.S.S. (black circles, taken from I Aharonian et alj d2007l) ). as 
well as the EBL corrected, intrinsic source spectrum (grey open 
diamonds) with a hard photon index, which implies an inverse 
Compton emission peaking at very high frequency (> 10 27 Hz), 
is shown. The solid line represents an SSC model that can de- 
scribe the intrinsic synchrotron and inverse Compton emission 
of 1ES 0229+200, and the dashed line represents the absorption 
by the EBL (details about these models are described in the text). 



0229+200 was also studied in the Bessel U, B, and V-bands wit h 
the Nordic Optical Telescope (NOT) bv lHvvonen et all (120071) . 
The results are mh os t,B = 18.75 mag, mh os t,u = 18.83 mag, 
and mhost.v = 17.58 mag with half light radii r e B = 5.65", 
r e u = 2.75", and r e v = 4.90". In order to correct for the host 
galaxy light, a de Vaucouleurs profile of the galaxy was assumed 
and the measured brightness w as transforme d to that of an aper- 
ture of 4" using equation (4) oflYoung ( 1976). ATOM photome- 
try was performed with a 4" aperture. The resulting host-galaxy 
corrected fluxes would result in an unphysical bump in the V 
band in the SED (in vF v ) since the calculated influence of the 
host galaxy in the V band is only 30%, while in the R and B band 
it is 90% and 30%, respectively. Hence, the influence of the host 
galaxy in the R, B, V, and U filters were also cal culated using a 
spectr al template for a nearby elliptical galaxy by Fukugit a et alj 
(1 19951) . Her e, the R-band magnitudes an d the half-light radius 
detected by iFalomo & Kotilainenl (1 19991) were used. The influ- 
ence of the host galaxy is then ~ 90%, ~ 80%, and ~ 57% for R, 
V, and B, respectively, which are the values that we used to cor- 
rect the measured fluxes (shown in Fig.[4]i. The host-galaxy cor- 
rected flux in the R-band is compatible wi th th e detected nucleus 
magnitude of lFalomo & Kotilainenl (1 19991) and lUrrv et all (2000) 
as expected from the absence of variability. The host galaxy in- 
fluence in the UVW1 and UVM2 bands is unknown. Figure [4] 
therefore shows two values, connected by a bar. The upper ones 
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correspond to the measured values corrected only for extinction, 
the lower ones also assume a correction for the host galaxy of 
30% (the value derived in the adjacent U band). 

In an independent check, the spectral slope was extracted 
using the Sloan Digital Sky Survey (SDSS) observations (five 
bands taken simultaneously). The resulting slope is identical to 
the one measured by ATOM. Since the data do not match the 
epoch of the ATOM observations, the absolute fluxes were not 
considered. 

5.2. Galactic extinction/absorption correction 

Since the host galaxy magnitudes were not extinction corrected, 
the influence of the host galaxy was first subtracted. The extinc- 
tion correction was thereafter applied to the AGN light. 

The measured UV fluxes w ere corrected for du st absorp- 
tion using E(B-V)=0.135 m ag dSchlegel et alj [l998h and the 
A,i/E(B — V) ratios given in Seaton ( 19791) . resulting in a cor- 
rection of 70%, 57%, 48% for the UVM2, UVW1 and U-band, 
respectively. For the optical filters, the values for the extinc- 
tion wer e derived from th e interstellar reddening curve and table 
given bv lZombeckl (fl 990FI The correction of Nn absorption was 
applied as described in Section 3. 

5.3. Historical multi-wavelength data 

GALEX observed 1ES 0229+200 at October 29, 2007 in the 
far UV (152.8 nm) and near UV (227.1 nm) leading to 
a measured flux of 31.68 /iJy and 50.36 yuJy, respectively, 
taken from GalexView^. The measured near and far UV fluxes 
were corrected for dust absorption using E(B-V)=0.135 mag 
dSchlegel etal.l 19981) and the A A /E(B - V) ratios given in lSeatonl 
d 19791) . A correction of 64% and 68% resulted for the far and 
near UV bands, respectively. Owing to a lack of information 
about the host galaxy influence in these wavebands (this influ- 
ence should be very small compared to the extinction correc- 
tion), the measured fluxes are shown as upper limits for the in- 
trinsic synchrotron spectrum of 1ES 0229+200 (shown as grey 
arrows in Fig. |4j». 

In 2006, 2008, and 2009 Integral also observed 1ES 
0229+200 several times for approximately 66 ks with ISGRI 
(17 - 80 keV) and 61 ks with JEM-X (3 - 10 keV), but the 
source appears faint and too few photons were detected so that 
no reasonable light curve or spectrum could be extractecQ 

In historical snaps hot observations of VL A at 6 cm in 1992 
fluxes of 41.5 mJy (ISchachter et all 1 19931) and 49.09 mJy 
dPerlman et al.ll 199a) were detected (shown as grey open circles 
in Fig . [5j> ■ With 6 cm VLBA observations th is flux was resolve d 
into a core of 22.7 mJy and a jet of 7.7 mJy dRector et al.l l2003). 
In the same study, the core flux at 1 .4 GHz was detected by VLA 
observations to be 51.8 mJy, and 1ES 0229+200 was found to 
have curved jets to the north and south with extensions of 30". 
The core fluxes are shown as grey filled circles in Fig. [5] 

5.4. SSC model 

Katarzvhs ki et alJ ((2006) demonstrated that synchrotron spectra 
with a high minimum Lorentz factor can explain a very hard 
inverse Compton spectra. The simultaneously obtained data in 
the optical, UV, and X-ray bands cover the whole synchrotron 

8 online version http://ads.harvard.edu/books/hsaa/toc.html 

9 http://galex.stsci.edu/GR6 

10 see e.g. http://www.isdc.unige.ch/heavens_webapp/integral 



emission. This allows an empirical determination of the mini- 
mum Lorentz factor in 1ES 0229+200. The high accuracy of the 
XMM-Newton measurements constrains the synchrotron spec- 
trum very efficiently. 

The photon index F = a + 1 of the X-ray spectrum gives a 
direct estimate of the spectral index n of the electron distribution 
N(E)dE oc E~"dE with the relation n = -1 -2a. A broken power 
law with a canonical cooling break of An = 1 is assumed for the 
electron distribution. The variability of 1ES 0229+200 detected 
in the X-ray emission places constraints on the maximum size 
of the emission volume considered in the SSC model which is 
inferred to be R < (Af/sec) xDxcr;Dx2.6x 10 14 m. 

The peaks in the SED a re commonly explained by lep- 
tonic models, such as those of iMarscher & Gear! d 19851) . as syn- 
chrotron and inverse Compton (IC) emission from a popula- 
tion of relativistic electrons up-scattering their self-produced 
synchrotron ph otons (synchrotron self-C ompton models (SSC)). 
The code of iKrawczvnski et all d2004l) for a one-zone SSC 
model was used to describe the intrinsic emission of 1ES 
0229+200. For this model, a spherical emission volume of ra- 
dius R, moving with a bulk Doppler factor D towards the ob- 
server and a magnetic field B was assumed. The electron distri- 
bution was described by a broken power-law between minimum 
and maximum energy. 

The X-ray spectrum can be fitted well by a single power- 
law model over 1.5 orders of energy. The extrapolation fits 
the absorption and host-galaxy corrected UV data, suggesting 
that there is a single power law slope for the energy range 
0.005 - 100 keV. The absorption and host galaxy corrected UV- 
optical emission below this range is significantly steeper and 
strongly constrains the minimum energy of the electron distri- 
bution function, yielding a very high value. The well-determined 
spectral index in the X-ray regime has to describe the low-energy 
spectral index of the electron distribution n\ = -2.6. The un- 
cooled electron spectral slope hence has «2 = -3.6 follow- 
ing the canonical break. The radius was kept as the maximum 
value obtained from the variability timescale (R = 1 X 10 16 m). 
This radius is an upper limit and is fixed independently of the 
SED modelling. Smaller values o f radii would also be c onsis- 
tent with the variability constraint. iTavecchio et al. (2009) used 
a smaller value in their attempts to describe the SED compiled 
in their study. The SED shown in Fig. 5 cannot be reproduced 
with a radius (R < 1 x 10 15 m) (see below). To account for 
the high energy peak of the IC emission, the underlying elec- 
tron distribution must be very narrow. Together these constraints 
are best met with the model parameters E m [ n = 2 x 10 11 eV 

(ymin = 3.9 X 10 5 ), £ brea k = 3.2 X 10 13 eV (ybreak = 6.2 X 10 7 ), 

Zi max . 1 x 10 14 eV (y max = 1.9 x 10 8 ), and a magnetic 
field of B - 3.2 x 10 5 G. The Doppler factor of D=40 was 
chosen as the smallest one possible to reproduce the spectra. 
This model can describe the simultaneously obtained spectra 
in the synchrotron regime, the long-term hard X-ray spectrum 
by Swift/BAT, and the non-simultaneous, intrinsic TeV spec- 
trum. The peak frequency of the synchrotron and IC emission 
based on the described model are v sy pea k = 3.5 x 10 19 Hz and 
vic.peak = 1.5 X 10 27 Hz. The peak fluxes show a slight Compton 
dominance with v S y,peakf v sy , peat = 9.7 x 10~ 12 erg cirT 2 s _1 and 

vic.peak-Fvic pe* = !•! x 10~ u erg cirT 2 s _1 . This model is shown 
in Fig. [5] as a solid line, while the dashed line represents the 
absor ption by the EBL using the model of iFranceschini et"aT1 
(2008). This model is in good agreement w ith the EBL absorp- 
tion model used in lAharonian et al.l (|2007) to correct the ob- 
served TeV spectrum (shown in Fig. [5j. 
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The cooling break is at very high energies such that a single 
power-law model for the electron distribution with n = -2.6 fits 
the data equally well and is an alternative model to describe the 
measured spectra of the SED. 

An earlier a ttemp t to utilize the suggestion of 
IKatarzvnski et al.l d2006l) was reported by pTavecchi o et al.l 
(2009)! They assumed a narrow single power-law electron 
distribution with a high minimum Lorentz factor to fit the data. 
As in the modelling presented here, the slope of the synchrotron 
emission below y m i n is given by F sy ~ v 1/3 , leading to a slope 
for the IC part that is similar to that is required to fit the 
TeV spectrum. Attempting to fit the overestimated X-ray flux, 
iTavecchio et al.l 12009) derived a different set of parameters, 
including a higher (factor 10) magnetic field and a smaller 
radius. The initial model presented above used an upper limit 
to the radius, which was not optimized in the SED fitting. In 
an attempt to explore models with smaller radii, attempts were 
made to reproduce the SED with smaller radii. The precise 
determination of the low energy cutoff in the UV range and the 
well-determined spectral shape in the X-ray range can only be 
reproduced by invoking a very high Doppler factor (D> 100). 
While providing an adequate fit to the synchrotron range, they 
do not reproduce the inverse Compton emission peak well. 

The cut-off of the synchrotron emission above 100 keV rep- 
resents the maximum energy of the electron distribution. No 
change in spectral shape could be detected in the 5w//f/BAT 
spectrum. Therefore only a lower limit to the maximum cut-off 
energy could be determined empirically. The precise measure- 
ment of the maximum energy would require a soft y-ray de- 
tection. In the IC part of the spectrum, the maximum energy 
cannot be determined emprically either because any increase in 
£max does not result in a change in the IC spectrum because of 
Klein-Nishina suppression. Assuming that the high energy cut- 
off of the synchrotron emission is at the high energy end of the 
Swift/BAT spectrum (~ 100 keV), this results in an electron dis- 
tribution that extends only over five orders of energy. 

Changes in the peak fluxes of the synchrotron emission (such 
as the difference between the XMM-Newton and the Swift/XKY 
spectrum, as shown in Fig. [4]i were also tested. Any changes in 
the parameters of the emission volume corresponding to an in- 
crease in flux by a factor of ~ 2 yield a change in the inverse 
Compton peak flux of a factor of ~ 2 as well. This implies that 
changes in the VHE fluxes in the TeV energy range are also ex- 
pected to occur in 1ES 0229+200. 

Alternative models were also explored. Apart from the at- 
tempts to describe the SED by SSC emission from a very com- 
pact source, other assumption would also lead to models involv- 
ing a high Doppler factor. Assuming that the derived curvature in 
the optical regime would be represented by the break of the elec- 
tron distribution with n\ = -1.6, «2 = -2.6, a separate model 
can be constructed. This scenario cannot be excluded by our si- 
multaneously measured SED above 10 14 Hz, but would require 
a large Doppler factor of about D=100 in order to account for 
the high energy peak of the IC emission. This alternative model 
also underpredicts the radio flux of the resolved core by a large 
factor, e.g. ar ound 70%. Interesti ngly, the radio flux of the jet, as 
measured by iRector et al.l d2003l) . could be represented by such 
model, indicating that the whole synchrotron emission originates 
from the resolved jet. Another possible way of explaining the 
high IC peak using such an electron distribution function, but 
with a lower Doppler factor, would be to invoke an additional 
external Compton contributi on resulting from t he scattering of 
the CMB as discussed in iBottcher et all d2008l) for 1ES 1101- 
232, a distant TeV blazar with a hard TeV spectrum. This model 



would decouple the X-ray and TeV emission because most of the 
TeV emission would result from the interaction with the CMB 
further away from the synchrotron emission region of the jet. 
This model would imply no VHE variability on short timescales. 
The effect of the CMB scattering depends on the distance to the 
source, so that it is more likely to be detectable in distant sources 
such as 1ES 0229+200 and 1ES 1101-232. 



6. Conclusions 

The X-ray spectral information that have been presented here 
cover all available X-ray data up to October 2010. Together with 
the simultaneous UV and optical observations, they yield a very 
good coverage of the synchrotron emission. Our detailed study 
of the high quality XMM-Newton spectrum has inferred an ab- 
sorption higher than the nominal Galactic column density, which 
could be intrinsic to the source or caused by Galactic excess ab- 
sorption. 

The host-galaxy, extinction-corrected optical and UV fluxes 
have been shown to provide strong evidence that the cut-off of 
the low energy part of the synchrotron emission is located be- 
tween the optical and X-ray regime. Therefore, the minimum 
energy of the electron distribution has to be rather high to ac - 
count for this cut-off. As suggested by Katarzvhs ki et al.l (|2006), 
an electron distribution with a high minimum Lorentz factor 
is needed to reproduce a hard TeV spectrum. A narrow elec- 
tron distribution, as indicated by the high minimum energy, re- 
sults in a hard intrinsic VHE y-ray spectrum as deduced for 
1ES 0229+200. This hard VHE spectrum also shows that the in- 
verse Compton peak is at very high energies (> 10 27 Hz), hence 
the initial electrons are very energetic. 

An earlier a ttemp t to utilize the suggestion by 
IKatarzvnski et al.l (120061) was reported by ITavecchio et al.l 
(120091) . assuming a high minimum Lorentz factor. The main 
differences between the model presented here and their attempts 
are caused by the new, precisely determined low-energy cut-off 
in the optical regime and the simultaneously determined X-ray 
spectrum. 

A broken power-law could be an alternative description to 
the additional absorption found in the X-ray spectra. However, 
the low energy extrapolation of the X-ray spectrum would not fit 
the optical range and is therefore disfavoured. 

The hard X-ray spectrum up to 15 keV together with the 
long-term spectrum by ^w/ff/BAT have been found to show that 
the synchrotron peak is extended up to ~ 100 keV without any 
significant cut-off in the X-ray spectrum. 

1ES 0229+200 is defined as a high-frequency peaked BL 
Lac object, and the measured synchrotron emission peaks at 
higher frequencies (> 100 keV) than usual for HBL and be- 
longs therefore to the class of extreme blazars. The exact peak 
frequency cannot be determined, since the hard X-ray spectrum 
does not show any significant cut-off. Assuming that the high en- 
ergy cut-off of the synchrotron emission is at the high energy end 
of the Swift/BAil spectrum (~ 100 keV), the underlying electron 
distribution extend only over five orders of magnitude in energy, 
which is a narrow range. For no other extreme blazar has this 
range been demonstrated to be so narrow. 

The low energy part of the synchrotron emission has not been 
found to show significant variation over four years in the optical 
R band. Only minor variation has been detected in the ultraviolet 
flux during 2008 and 2009. The 2-10 keV X-ray flux instead 
varied by a factor of « 2 within around 20 days in 2009. In 
the high energy part of the synchrotron emission, measured by 
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SwiftfB AT, variations in this amplitude and timescale could not 
be detected because of the longer integration time needed. 

Several blazars are known as so-called extreme blazars 
with synchrotron emi ssion extending to high energies 
dCostamante et all 1200 lb . It has been found that Mkn 501 
revealed a very high energy peak of synchrotron emission 
around 100 keV in a flare with a detected larg e shift in the peak 
frequency compared to previous observations dPian et alj fl998). 
The monitoring by BeppoSAX of 1ES 2344+514 has identified 
huge changes in the synchrotr on peak frequency w ithin one year 
with different spectral shapes dGiommi et al.ll2000t) . INTEGRAL 
observations of 1ES 1426 +428 display a sy nchrotron peak 
frequency around 100 keV (Wolt er et alJ l2008). These extreme 
blazars have similar synchrotron peak frequencies to 1ES 
0229+200, but their very high energy spectra in the TeV range 
are much softer than for 1ES 0229+200. Hence, 1ES 0229+200 
has the highest IC peak frequency among the extreme blazars 
known up to date. 
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Table 3. X-ray sources detected in the XMM-Newton observation on 1E S 0229+200. The ir possible counterparts in the ra dio, IR, 
and optical regimes were searched in the Guide Star Catalog (GSC) dLasker et alJl2008h . the IRAS faint source catalog (|Moshir 
[l990i and the NRAO VLA Sky Survey (NVSS) dCondon et alj[l998h . 



